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Abstract 

We describe how exactly the intra-cavity helds in a dual recycling cavity build up 
their power before achieving a steady state value. We restricted our analysis here to 
interferometers with lossless mirrors and beam-splitter. The complete series representation 
of the intra-cavity lights at any stage of evolution in non-steady state have been presented. 
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Laser interferometric gravitational wave observatories [1] are currently being designed 
or constructed in several countries (LIGO, VIRGO, GEO, AIGO). The interferometer, in 
essence, is a rather elaborate transducer from optical path difference to output power. 
So, monitoring the change in output power, it would be possible to detect the changing 
curvature of spacetime induced by the passage of a gravitational wave. In order to reduce 
the shot noise level, an integral feature of these detectors will be the use of high power laser, 
in conjunction with variants of optical technique known as light recycling. The first one 
of these techniques is called power recycling[2] , in which, at a dark fringe operation of the 
interferometer, the outgoing laser light is recycled back into the interferometer by putting 
a mirror in front of the source, thus enhancing the laser power. When a gravitational 
wave passes through an interferometer, it modulates the phase of the laser light, thus 
producing sidebands which travel towards the photodetector [3]. So, another variant of 
optical technique called dual recycling incorporates a second recycling mirror placed in 
front of the photodetector as shown in Fig.l. This arrangement can store signal sidebands 
for suflBciently long time to allow optimum photon noise sensitivity within a restricted 
bandwidth [4]. 

In this note, we describe, how with the help of Mathematica[5] , one can do calculation 
for the evolution of intra-cavity fields in the two recycling cavities of a dual recycled 
interferometer and finally achieve the expressions for the intra-cavity fields at the steady 
state condition. This may find its application in investigating various kinds of transient 
states in such interferometers. 

We first introduce the notations used here. In all the following figures, rays of light 
labelled by a, 6, x, y represent the complex amplitude of the beam's electromagnetic field. 
Mirrors: In this note, we make the following assumptions on the mirrors and beam- 
splitter: (i) these are lossless, (ii) the substrate and the dielectric layers in mirrors are 
linear media, (iii) mirrors are time-reversal invariant, (iv) any mirror has a reflection 
symmetry with respect to exchange of inputs and outputs on the same side of the mirror, 
(v) there is no differential time delay between the transmitted and reflected lights and also 
the overall time delay across the mirror is negligible. One can then arrive at the following 
simple input-output mirror relation (refer to Fig. 2): 



where t and r represent the transmission and reflection coefficients of the mirror. For a 




(1) 
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lossless mirror, | t p + | r p= 1. We consider all the beam-splitters to be of 50:50 type. 
Fabry-Perot Cavity: The flowchart of the light paths in such a cavity has been shown 
in Fig. 3. Light beams ao and bo enter the cavity from corner mirror (CM) and end mirror 
(EM) sides respectively. The beam x represents the intra-cavity field that grows up in 
every bounce it gets from CM or EM, provided it adds up in phase with the incoming light 
there. The index j represents the number of round-trips of the x-field inside the cavity. 
The circle with the quantity L inside it represents the change in phase of the intra-cavity 
light due to the traversal of cavity of length L. After a sufficient number of bounces, 
the intra-cavity field, x, achieves a steady state when it stops growing any more (i.e., the 
change after any bounce would be negligibly small) . 

Dual Recycling Cavity: The corresponding fiowchart for a single delay-line dual- 
recycled cavity has been shown in Fig. 4. To understand this figure one may also refer 
to Fig.l. The circle with the quantity li or I2 inside it represents the change in phase of 
the intra-cavity light due to traversal of power recycling cavity (PRC) of length li between 
the beam-splitter (BS) and the power recycling mirror (PRM) or Signal Recycling Cavity 
(SRC) of length I2 between BS and signal recycling mirror (SRM) respectively. The elec- 
tromagnetic field ao represents the coherent beam of laser light entering through PRM. 
The field bo is in vacuum state in ordinary interferometers, but as suggested by Caves [6], 
in advanced interferometers, this may be made to be in squeezed vacuum state to reduce 
the shot noise. The fields X4 and represent intra-cavity light beams in the Power Re- 
cycling Cavity (PRC) and Signal Recycling Cavity (SRC) respectively. We represent the 
amplitude transmission and reflection coefficients of the PRM (SRM) by ^1(^2) and ri(r2) 
respectively. The box with the label 'Propagation in Arms' represents the matrix relation: 



where ko = ujq/c and Uq, La, c are the circular frequency of light beam, the arm length 
of the interferometer and the velocity of light respectively. The quantity 9 represents the 
differential phase shift of light in the two arms of the interferometer. The end mirrors of 
dual recycling cavity are assumed to be perfectly reflecting ones. 

A small programme in Mathematica can be easily written to implement this recursive 
procedure. After a sufficient number of bounces, j (say, 10), one can look at the expressions 
for X4[j] and y4[j]. We have set a;o[0] = yo[0] = in the final expressions of X4[j] and y4[j] 
after finishing the recursive calculation because the initial intra-cavity fields were zero. 
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Now, grouping the related terms together, we can guess the formation of a series in each 
case. 

The series obtained by us are as foUows: 
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where 



Qo = ex.p[iko{2La + h + h)] 
Qi =exp[iko{2La + 2h)] 
Q2 ^exp[iko{2La + 2l2)] 



(3c) 



are the phase factors which satisfy Qq = QiQ2- The upper hmits of the integers k, n and 
m are represented by K, N, and M respectively. For any value of j, we obtain all the 
terms in the series that are characterised by sets of values of these limits satisfying the 
following relationship: 

j = 2K + M + N. {3d) 

So, in the complete expression for X4, say, there would be several terms each involving the 
product of three summations over k, m and n, one term for each combination of M 
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and N satisfying Eq.(3d). So, as j increases, the number of triple products of summations 
in Eqs.(3a) and (3b) also increase. The explanation of how these upper limits arise is 
provided below. 

In this connection, interested readers may also like to know the physical interpretation 
of all the terms appearing in series (3a) and (3b). One may note that any light gets 
multiplied by either (z sin 9) or (cos 9) depending on whether it changes its cavity (from 
PRC to SRC and vice versa) or not after passing through the BS and one arm. 

For example, we may consider the series in the second term of the coefficient of ao in 
the intra-cavity field x^lj]. Any term in this series represent the contribution of ao light 
in X4 after it completed {k + n) number of round trips in PRC and {k + m) number of 
round trips in SRC. Since j represents the sum of the maximum number of round trips 
in two cavities, this explains why we obtain the relation (3d) for the series described in 
Eqs.(3a) and (3b). The number of possible ways of completing such a trip would be 
^n+k(j^^^y^m+k-iQ^y_ As an easy example, let us take the contribution corresponding to 
A; = 0, n = 2, m = 1, i.e., — 2rir2ti cos^sin^ 6, which may be expanded as 



This means that this term has been contributed by two parts of ao light which made two 
round trips of PRC and one of SRC in the following two ways: 

(i) ^PRM; transmitted inside ->BS, arm ^SRC ^SRM; reflected ^BS, arms ^PRC 
^PRM; reflected ^BS, arm ^PRC. 

(ii) ^PRM; transmitted inside ^BS, arm ^PRC ^PRM; reflected ^BS, arms ^SRC 
^SRM; reflected ^BS, arm ^PRC. 

The intra-cavity fields, and achieve their steady state value as j tends to infinity. 
Now, identifying 



etc. and after doing some algebraic manipulations, we finally arrive at the steady state 
expressions: 



(ri cos 9) {i sin 9)r2 {i sin 9)ti + {i sin 9)r2 (^ sin 9) (ri cos 9)ti . 
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where 

X = 1 - (riQi + r2Q2) cos 9 + rir2QiQ2- (5c) 

The steady-state expressions (5) can, of course, be calculated very easily by assuming that, 
under steady-state condition, the change in the intra-cavity fields after every round-trip of 
the cavities is negligibly small [7]. One can then write the following two equations which 
can be solved to obtain the above-mentioned expressions. 

a;4[steady] =aQtiQi cos 9 -\- 0:4 [steady] riQi cos 9 

+ ibot2Qo sin 9 + zy4 [steady] r2Qo sin 9, (6a) 
1/4 [steady] —bot2Q2 cos 6 + ^4 [steady] r2Q2 cos 9 

+ iaotiQo sin 9 + ia;4 [steady] riQo sin 9. (66) 

The intra-cavity fields can achieve their maximum values at the resonance condition of the 
dual recycling cavity, i.e., when li = I2 and Qq = Qi = Q2 = 1- This implies that the 
intra-cavity fields X4 and ^4 add up in phase with the incoming light. For a sufficiently 
low value of r2, this condition also represents the broadband mode of operation of a dual 
recycled cavity [4] . The expressions for the case of power recycling (when no SRM is used) 
can be easily obtained by setting r2 = 0. Similar exercise can of course be carried out for 
interferometers incorporating Fabry-Perot cavities in the arms. 

For exactly equal armlengths of the interferometer (^ = 0), the dual recycling cavity, 
at any frequency, gets decoupled into two equivalent three mirror cavities. If Qi = 1, the 
laser light, ao becomes resonant with the three mirror power recycling cavity (PRM and 
two EMs), whereas the sidebands of laser light (generated by the gravitational wave [3]) 
travel down the SRC. One can make one of these sidebands (say, ujq + ujg) resonant with 
the three mirror signal recycling cavity (SRM and two EMs) by adjusting I2, such that 
Q2ex.p[+iujg{2La + 2l2)] = 1. This is called the narrowband mode of operation. Obviously, 
under this mode of operation, neither the bo light nor its sidebands ( which, after getting 
created, travel down the PRC) can become resonant with any of the cavities. 

Gravitational radiation from a coalescing compact binary system has a chirp waveform 
which continually increases in frequency and amplitude as the two stars spiral in towards 
each other. To observe this broadband waveform, dynamical dual recycling technique[8] 
was proposed, which incorporated two complementary techniques to improve the sensitivity 
of observation: (i) narrow band observation to achieve high signal to noise ratio at a 
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given frequency, and (ii) tracking the chirp frequency by adjusting the tuning frequency 
by changing I2 as long as the signal lies within the detector bandwidth. To implement 
this technique, it is important to understand how exactly the signal (i.e., the sideband 
that is tracked down) builds up its power in the presence of errors in tuning the SRC. A 
programme simulating this technique in the way described above, therefore, may come to 
be handy and useful. 

In actual LIGO type interferometers, the Fabry Perot cavities in the arms may be 
physically different in length by several meters, and the storage times in the arms different 
by as much as 1.0% due to differing mirror transmissions or deviations from the operating 
point [9] . It may be noted that this way of doing the calculation exactly tracks down various 
light beams and thus provides an answer to the question 'how many times a particular 
beam traverses through any particular arm or cavity before a specific point of time in 
a transient state'. The arbitrary phase noise (mainly caused by various nonlinear effects, 
e.g. scattering, thermoelastic effects on the mirrors etc.) induced on the beams in different 
arms may lead to leakage of laser light to the SRC. Also, an important issue is the control 
of the resonance conditions in the various cavities in these type of interferometers made 
of suspended mirrors. For this purpose it is useful to compute the transient behaviour 
of the light amplitudes at various locations in the interferometer for given motions of 
mirrors. This particular algorithm will help in investigating such kinds of transient state 
of the interferometer and in tracking down the contribution of different cavities in the total 
phase noise in an effective way. Investigations are currently in progress in this line and 
results will be communicated in future. 

One point of caution may be worth noting here, although that is not important for 
the present designs of the interferometric gravitational wave detectors. If the Fabry Perot 
cavities in the arms of the interferometer are significantly different in length, then within 
a given time period, the number of bounces of light in one cavity will become different 
from the number of bounces in the other. The analysis in the form presented here will 
break down in that case because the iteration parameter j will differ for different beams 
depending on the path history. 

We point out again that the series compaction to arrive at Eqs.(3a) and (3b) was per- 
formed by manual inspection of the recursive series expansion produced by Mathematica. 
Employing Mathematica to find these automatically is, in general, impossible since there 
may be infinite number of functions with identical series expansions up to a certain degree. 
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Abbott [10] has studied different aspects of this Mathematica problem and showed a method 
of compacting series using generating functions, if the corresponding recurrence relation 
is available. Investigation is in order to find a way to combine the complicated recursive 
process of the present problem with the method described by Abott[10] for simpler cases, 
with the purpose of employing Mathematica to arrive at Eqs(3a) and (3b) automatically. 
It should be noted, however, that this is not a problem for the applicability of the method 
of calculation described in this paper to realistic situations where one has to do numerical 
computation. 

It is a pleasure to thank S.V. Dhurandhar and B.S. Sathyaprakash for discussions and 
Paul Abbott for taking active interest in the problem of 'series compaction' on my request. 
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Figure Caption 

FIG.l : Optical arrangement for dual recycling. (BS) - Beam Splitter, (EM) - End Mirror, 
(PRM) - Power Recycling Mirror, (SRM) - Signal Recycling Mirror, (PD) - Photode- 
tector. The light beams: ao - input laser beam, bo - ordinary or squeezed vacuum 
state (see teod); xi, x/^- intracavity fields in power recycling cavity; j/i, y4 - intracavity 
fields in signal recycling cavity; cci, X2i 1/2, J/3 - light fields inside the arms; ai, bi - 
output beams through PRM and SRM respectively. 

FIG. 2 : Mirror and light beams, gq, bo - input beams; ai, 61 - output beams. 

FIG.3 : Flowchart of a Fabry-Perot cavity at non-steady state. CM - Corner Mirror, EM 
- End Mirror. The circle with L inside represents phase change due to traversal of 
cavity length, L. 

FIG. 4 : Flowchart of a single delay-line dual recycled cavity at non-steady state. PRM - 
Power Recycling Cavity, SRM - Signal Recycling Cavity, BS - Beam Splitter. Circle 
with li or I2 inside represents phase changes due to traversal of length, li or I2 between 
BS and PRM or SRM respectively. 
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